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Abstract 

A simple model of cosmic ray electron acceleration at the jet boundary (Os- 
trowski 2000) yields a power-law particle energy distribution of ultrarelativistic 
electrons with an energy cut-off growing with time and finally, a growing particle 
bump at the energy, where energy gains equal radiation losses. For such elec- 
tron distribution, in tens-of-kpc scale jets, we derived the observed time varying 
spectra of synchrotron and inverse-Compton radiation, including comptoniza- 
tion of synchrotron and cosmic microwave background photons. Slowly varying 
spectral index along the jet in the 'low- frequency' spectral range is a natural 
consequence of boundary layer acceleration. Variations of the high energy bump 
of the electron distribution can give rise to anomalous behaviour in X-ray band 
in comparison to the lower frequencies. 
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1 Introduction 



Particle acceleration in a velocity shear layer was proposed in the early eighties by 
Berezhko with collaborators and then gradually developed (cf. a review by Berezhko 
1990). Action of such a mechanism at relativistic shear layers occurring at side bound- 
aries of relativistic jets was discussed by Ostrowski (1990, 1998, 2000). The mechanism 
was considered to provide ultra high energy cosmic rays and to be an important factor 
influencing the dynamics of relativistic jets in extragalactic radio sources thereby form- 
ing cosmic rays cocoons. Ostrowski & Sikora (2001) proposed that such acceleration 
mechanism can provide a cosmic ray proton population being a substantial pressure 
factor in FR II radio source lobes. 

The existence of the boundary layer with a velocity shear was first suggested by the 
observations of radio and optical jets in FR I sources. Owen et al. (1989) concluded that 
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the jet morphology in M87 is dominated by the boundary layer more than by the shock 
structure. Constancy of radio-to-optical spectral index along the jet and similarity of 
radio and optical maps with emissivity peaking near the jet's surface strongly suggested 
acceleration of radiating electrons in situ at the boundary layer. Later radio and optical 
polarimetry of M87 (Perlman et al. 1999) confirmed the complex structure of its jet 
containing a shear layer with a highly ordered magnetic field parallel to the jets axis. 
A map of the radio spectral index of another well known BL Lac, Mkn 501, revealed a 
boundary region of its jet with an unexpectedly flat spectrum (Edwards et al. 2000), 
suggesting action of some reacceleration mechanism at the jet surface. 

Some FR 11 radio sources also reveal spine - shear layer jet morphology. For in- 
stance radio observations of 3C353 (Swain et al. 1998) was interpreted in terms of a 
Doppler-hidden relativistic spine surrounded by the slower moving boundary region. 
VLBI observations of another FR II object 1055-^018 (Attridge et al. 1999) revealed 
a fragmentary but distinct boundary layer with the mean longitudinal magnetic field 
present where interactions with surrounding matter were strongest (i.e., where the jet 
bended). Optical and radio observations of the quasar 3C 273 (Jester et al. 2001) 
showed smooth changes of the radio-to-optical spectral index along the jet. This 
strongly suggested reacceleration of the radiating electrons taking place not in the 
specific shock-like regions, but within the whole iet. 

One should note that our knowledge of physical conditions at the relativistic jet 
boundary is still insufficient. Cosmic ray viscosity and radiation viscosity in a boundary 
layer with a velocity shear influencing hydrodynamics of the jet and particle distribu- 
tion along its surface was studied by Earl et al. (1988), Jokipii et al. (1989) and 
Arav & Begelman (1992), respectively. On the other hand, 3D simulations revealed 
the presence of a jet shear layer with a high specific internal energy and highly tur- 
bulent, subsonic and thin cocoon surrounding the relativistic fiow (Aloy et al. 1999). 
Simulations of radio emission from such jets (Aloy et al. 2000) show radial structures 
in both intensity and polarization. 

2 Acceleration of Cosmic Ray Electrons 

Let us assume that a thickness of the transition layer between the jet and the surround- 
ing medium, D, is comparable to the jet radius Rj ^ 1 kpc (cf. Owen et al. 1989, Swain 
et al. 1998), and that the magnetic field B frozen into tenuous plasma at and outside 
the boundary is on average parallel to the flow velocity U. Acceleration of energetic 
particles injected into the boundary layer results from its scattering on irregularities 
of the magnetic field in a turbulent medium with a perpendicular to U mean velocity 
gradient. In the considered case of the mean magnetic field aligned along the jet in 
a highly turbulent boundary layer plasma, a particle gyroradius can be considered 
instead of a particle mean free path A, and cross- field radial diffusion provides a particle 
escape mechanism from the acceleration region. For particles with energies satisfying 
Tg <^ D, the acceleration within a finite thickness turbulent shear layer results from, 
both, the second order Fermi process in the turbulent medium and the cosmic ray 
viscosity. We consider the jet velocity to be comparable to the light velocity even on 
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kiloparsec scales, U ^ c, and a characteristic velocity of magnetic turbulence, V, to be 
comparable to the Alfven velocity for a subsonic turbulence and at least a few orders of 
magnitude smaller than U. The Fermi process dominates over the viscous acceleration 
when the particle energy is small enough to satisfy < D(y/U). Because of rapid 
energy radiation losses at high energies, the electrons are always expected to satisfy 
the required condition. Then, the acceleration time scale in a highly turbulent shear 
layer can be estimated (Ostrowski 2000) as 

T (? 

Tacc ~ (1) 

With anticipated parameters of the boundary layer - B ^ 10"^ G and V ^2 - 10® cm/s 
- the characteristic acceleration time scale (1) is T^cc ~ 10^ 7 [s], where 7 is the electron 
Lorenz factor. One should note, that equation (1) illustrates the most optimistic 
scenario for the acceleration process, with presence of non-hnear MHD turbulence at 
the scale of Vg. The time scale for particle escape from the acceleration region due to 
cross-field diffusion can be estimated as 

Tesc ~ , (2) 

where is the cross-field diffusion coefficient, k± = rjrg c/3, and 77 is a scaling numer- 
ical factor. The particle escape becomes important when the electron energy satisfies 
condition Tacc ^ ^esc ; i-G- whcu r-g > {?>/riY^'^ D iy / c). For the considered bound- 
ary layer parameters and a strong turbulence condition 77 ~ 1, equation (2) reads 
as Tesc ~ 6 • 10^^ 7~^ [s], and hence radiative (synchrotron) cooling with a time scale 
Tcooi ~ 8 • lO"*"® 7~^ [s] becomes important for the much lower electrons energies. There- 
fore, by comparing the acceleration time scale (1) with the time scale for radiative 
losses, one obtains the maximum electron energy, 7eg ~ 10®. Rapid radiative cooling of 
electrons with such high maximum energy is compensated by the continuous accelera- 
tion taking place within the whole boundary layer along the jet. For lower turbulence 
amplitude the acceleration process can proceed slower to yield a lower value for 7eg. 

In the present model, after switching on the injection of seed electrons with a small 
initial Lorentz factor 70, a time dependent spectral evolution of particles accelerated in 
the turbulent layer begins as a power-law distribution with a growing energy cut-off. 
Then it evolves into a characteristic shape consisting of two different components: a 
stationary flat power-law part at low energies, ricij) oc 7"^, finished with a nearly 
monoenergetic growing peak, or bump, at 7 = 7eg, due to an accelerated particles pile- 
up caused by losses (cf. Ostrowski 2000). Such electron spectrum can be approximately 
given by a formula using the Haeviside function, 0, and the Dirac delta function, 5. 
At energies above the injection energy, 70, the electron distribution 

ne(7, t) = 07"^ 6(7maa;(t) - 7) + Kt) ^(7 " leq) 0(7t " Icq), (3) 

where a and b{t) are normalization parameters (and b{t) is a growing function of time), 
7f is the highest energy the electrons can reach since the injection with radiation losses 
neglected, and -fmaxit) = min(7j, 7eq). 
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3 Radiation of electrons accelerated in a shear layer 



We consider two radiation processes taking place at the relativistic jet boundary: the 
synchrotron (SYN) emission of ultrarelativistic electrons with the energy spectrum (3), 
and their inverse Compton (IC) cooling. Synchrotron self-absorption can be neglected 
with the assumed optical depth for this process to be much less than unity for v > 
10^° Hz. Additionally, the present preliminary IC scattering calculations are limited to 
the Thomson regime assuming a step scattering cross section. A high energy spectral 
cut-off resulting from absorption of gamma rays due to photon-photon pair creation 
on a cosmic infrared background is expected to occur for distant sources (e.g. Renault 
et al. 2001). Therefore, the obtained below photon spectra above TcV energies cannot 
be directly compared to observations. In the case of IC emission of a tens-of-kpc scale 
jets, we study the synchrotron sclf-Compton (SSC) radiation and Compton scattering 
of external CMB photons (EC). With the assumed parameters, the magnetic field 
energy density at the jet boundary is -u^ ~ 4 ■ lO^"*^^ erg/cm^. The energy density of 
the CMB field in a source frame moving with a bulk Lorentz factor F at a redshift z 
is Ucmb = ■WcmfeT^ (1 + z)"^ , where an asterisk denotes quantities in our local rest frame, 
and u*.^j^i^ ^ 4 ■ 10~^^ erg/cm^. In the jet rest frame Ucmb can be comparable to ub, and 
both above processes can be significant even for small F and z. 

For the optically thin emission, radiation intensities of the power-law ('!') and the 
high energy bump ('2') electron components (c.f. Eq. 3) are approximatly 

J« (Z/, t) OC e{u,yn{t) - , t) OC e{u,yn{t) - l^) 6(7* - 7e,),(4) 

/£f )(^, t) OC u-"^ e(z/,,,(t) - u) , t) OC e(z.,,,(t) - U) 0(7* - 7e.), (5) 

I^Xv,t)<^v-^'^Q{v,,{t)-v) , 42)(z.,i)oci/2e(i/ec(i)-z^)e(7t-7e,), (6) 

where, for the assumed large scale jet parameters and neglecting the cosmological 
redshift correction, the maximum synchrotron, sclf-Compton and external Compton 
frequencies are Vsyn{t) ~ 10 7^„^(t)Hz, Vsscit) ^ min [v syn{t) , rnc^/h-fmax{t)] 7max(^) 
and I'ecit) ~ 10^^ F 7^„3,(t) Hz, respectively. In deriving equations (4) - (6) we used 
the standard theory of the synchrotron and the inverse- Compton radiation with simple 
approximations allowing evaluation of analytic expressions (e.g. Crusius & Schlick- 
eiser 1986 and Dermer 1995, respectively; for details see forthcoming paper Stawarz 
& Ostrowski 2001). The intensities (4) - (6) allow to describe the photon spectrum 
variations resulting from the particle distribution evolution given in equation (3). In 
the model, from the beginning of injection the accelerated electrons form the power-law 
spectrum with the energy cut-off growing with time. Then •jmax (^) = It and the result- 
ing emission is a power-law oc z/~^/^, with a growing cut-offs Ugynit), t^sscit) and Uec{t). 
When the maximum electron energy approaches 7eg, the high energy particle density 
peak forms and starts to grow. Then, an additional radiation components oc u^^^ (for 
the synchrotron radiation), oc and oc v^^^ (for the SSC emission), and oc v"^ (for 

the EC emission) appear and eventually dominate the spectrum at highest X-ray and 
7-ray frequencies. The spectral index « = — 1/3 results from the asymptotic expansion 
of the synchrotron emissivity function of the 'peak' electrons with the energy 7eq (c.f. 
Crusius & Schlickeiser 1986), and the index a = —2 corresponds to the Rayleigh-like 
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approximation of the anizotropic (in the relativistic fluid rest frame) CMB radiation 
comptonized by the monoenergetic electrons. 

Let us also briefly discuss the main implications of kinematic effects in the radiative 
jet output. For illustration, let us consider a long cylindrical jet of the radius Rj, 
moving with a constant flow Lorcntz factor Tj 10, surrounded by a boundary layer 
with a constant thickness D k, Rj ■^a Ikpc, (cf. Aloy et al. 2000), and the radial flow 
profile r(r). In the local source rest frame, the synchrotron, SSC and EC intensities 
are given by expressions (4), (5), and (6), respectively, but the observed intensities are 
It^v*, e*) = S'+'-{r, 9*)Uu*), and P^y^/^Ju*, 9*) = 5^+^{r, r) (z/*), where a 
is a spectral index given in the respective expression (4) - (6), and 6{r, 6*) is a Dopplcr 
factor for a sub-layer at a distance r from the jet axis and the jet inclination to the line 
of sight equal 9*. Because of a velocity shear, the boundary layer emission beaming 
pattern is different compared to the uniform beaming pattern of a spine radiation. For 
a given viewing angle 9*, boundary layer emission comes mostly from the region with 
r(r) = 1/ sin^*, where the Doppler factor is maximized. As a result, the observed jet- 
counterjet brightness asymmetry is decreased in comparison to the uniform jet with 
r = Fj , allowing the spine to be more relativistic in the large scale jets than traditional 
estimates suggest (cf. Komissarov 1990). 

The observed boundary layer emission, v* 9*), are plotted on Fig. 1 for four 

different viewing angles. Lines indicated by ti, t2 and correspond to different forms 
of the particle spectrum obtained in the considered acceleration model in successive 
times from the beginning of the acceleration process. For instance the line ti illustrates 
resultant sum of synchrotron, self-Compton and external Compton emission of the 
power-law electron component (index '1' in Eq. 4-6) with an energy cut-off '~fmax{t) = 
10^ < 7eq, while the line includes also the additional contribution from the high 
energy electron bump at '-/eq = 10^ (index '2' in Eq. 4-6). The presented spectra are 
integrated over the line of sight with the assumed linear dependence F = F(r) within 
the shear layer. 

4 Conclusions 

In the present paper, we study radiation of ultrarelativistic electrons accelerated at 
the boundary shear layer of the large scale relativistic jet. The considered acceleration 
process generate a characteristic particle distribution consisting of two components: 
the power-law section extended up to a more or less narrow pile-up bump, modeled 
here as the monoenergetic peak. A possible characteristic of such a cosmic ray popula- 
tion is its time variability, constrained by the not considered here, non-linear particle 
backreaction at the acceleration process. 

The acceleration within the jet boundary layer with inefficient particle escape due 
to radial diffusion could increas the energetic particle energy density Ug to extremely 
high values. However, when Ug becomes comparable to the ambient medium energy 
density, the accelerated particles can influence physical conditions within the boundary 
region. The involved electron pressure tends to expand the acceleration region rarefying 
the medium and, possibly, decreasing the efficiency of the acceleration process (cf. 
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Figure 1: The observed boundary layer emission, u* 9*), for four different viewing 

angles 9*, as measured near the jet, z = 0. Lines ti, t2, and correspond to 'jmaxi't) = lO^j 
lmax{t) = 10^, and -/max{t) = leq = 10^, respectively. 



Arav & Begelman 1992, Ostrowski 2000). Changes of the boundary layer density 
can result in changing several characteristic of the acceleration process, in particular 
shifting of the local value of and changing of the flow profile at the boundary. Such 
particle pressure backreaction, leading to local temporal variations of the high energy 
electron component, will generate fluctuations of the high energy radiation over the 
stationary power-law component emission at low photon energies. For the assumed 
physical conditions with '^^q = 10^, the process involves variations at X-ray and 7-ray 
bands with energies of ~ keV and ~ TeV, respectively. In our study of the tens-of-kpc 
scale jets, the involved time scales are ~ 10^° s, and, thus, can result in spatial kpc-scale 
brightness variations along the jet boundary. Much shorter scales are possibly involved 
in (sub) mili-arc-sccond jets, with much smaller (a factor of 10"^) spatial scales, much 
larger (a factor of 10^) magnetic fields, and stronger external radiation fields (study in 
progress now). 

As the described above radiative spectral components appear in a natural way in 

the discussed process, the stationary or time averaged spectra possibly formed due to 
the mentioned non-linear effects are also expected to posses such components. Then, 
one should observe X-ray and 7-ray excesses in electromagnetic spectrum of the jet. 
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as illustrated on Fig. 1. With the jet parameters considered above, such an X-ray 
flux above the extrapolated radio-to-optical continuum can explain the keV emission 
detected in several jets by Chandra. 

One should note, that our simple model, far from being complete, is supposed to 
describe main radiative properties of the boundary region under simple assumptions, 
and to provide restrictions on physical conditions within such layer. 
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